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ABSTRACT 

A  large-eddy  simulation  (LES)  model,  which  adopts  wave-averaged  equations  with  vortex  force,  is  used  to 
investigate  Langmuir  turbulence  and  ocean  boundary  layer  (OBL)  dynamics  in  high-wind  hurricane  condi¬ 
tions.  The  temporally  evolving  spatially  asymmetric  wind  and  wave  Stokes  drift  velocity  imposed  in  the  LES 
are  generated  by  a  spectral  wave  prediction  model  adapted  to  Hurricane  Frances  traveling  at  a  speed  of 
5.5  ms-1.  The  potency  of  Langmuir  turbulence  depends  on  the  turbulent  Langmuir  number,  the  wind-Stokes 
drift  alignment,  and  the  depth  scale  of  the  Stokes  profile  Ds  relative  to  the  OBL  depth  h.  At  the  time  of 
maximum  winds,  large-scale  vigorous  coherent  cells  develop  on  the  right-hand  side  of  the  storm  under  the 
inertially  rotating  winds;  the  Stokes  drift  velocity  is  well  tuned  to  the  surface  winds.  Much  weaker  cells 
develop  on  the  left-hand  side  of  the  storm,  partly  because  of  reduced  Stokes  production.  With  misaligned 
winds  and  waves  the  vertical  momentum  fluxes  can  be  counter  to  the  gradient  of  Stokes  drift,  and  the  cell 
orientation  tracks  the  direction  of  the  mean  Lagrangian  shear.  The  entrainment  flux  is  increased  by  20%  and 
the  sea  surface  temperature  is  0.25  K  cooler  on  the  right-hand  side  of  the  storm  in  the  presence  of  Langmuir 
turbulence.  Wave  effects  impact  entrainment  when  the  ratio  DJ\h\  >  0.75.  Because  of  wind-wave  asymmetry 
Langmuir  cells  add  quantitatively  to  the  left-right  asymmetry  already  understood  for  hurricanes  due  to 
resonance.  And  the  transient  evolution  of  the  OBL  cannot  be  understood  simply  in  terms  of  equilibrium 
snapshots. 


1.  Introduction 

The  interaction  between  winds,  waves,  and  currents 
ultimately  determines  the  fate  of  turbulent  mixing  in  the 
upper  ocean  boundary  layer  (OBL).  There  is  a  growing 
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body  of  large-eddy  simulation  (LES)  modeling  studies 
supported  by  observational  evidence  that  finds  surface 
waves,  and  in  particular  phase-averaged  wave-current 
interactions  and  stochastic  wave  breaking,  play  impor¬ 
tant  roles  in  the  dynamics  of  the  upper  ocean  boundary 
layer  (Sullivan  and  McWilliams  2010).  The  most  com¬ 
pelling  theoretical  model  of  wave-current  interactions  is 
embodied  in  wave-averaged  equations  and  their  math¬ 
ematical  representation  by  a  vortex  force  (Craik  and 
Leibovich  1976;  Leibovich  1983;  McWilliams  et  al.  1997; 
McWilliams  and  Restrepo  1999;  McWilliams  et  al.  2004). 
The  latter  is  found  to  generate  coherent  Langmuir  cells  on 
a  background  of  OBL  turbulence  that  enhance  turbulence 
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kinetic  energy  and  promote  entrainment  at  the  thermo- 
cline  (Skyllingstad  and  Denbo  1995;  McWilliams  et  al. 
1997).  Wave  breaking  can  modulate  the  formation  of 
these  coherent  structures,  and  breaking  energetics  in¬ 
crease  dissipation  rates  in  the  upper  OBL  by  several 
orders  of  magnitude  compared  to  classic  law-of-the-wall 
scaling  (Terray  et  al.  1996;  Sullivan  et  al.  2007;  McWilliams 
et  al.  2012).  The  regime  where  surface  waves  play  an  im¬ 
portant  role  in  the  dynamics  of  the  OBL  is  now  referred 
to  as  Langmuir  turbulence  (McWilliams  et  al.  1997). 
Langmuir  turbulence  is  created  during  the  passage  of 
storms,  cold  fronts,  and  tropical  cyclones.  These  inter¬ 
mittent  high-impact  weather  events,  which  help  define 
the  overall  ocean  structure  and  global  sea  surface  tem¬ 
perature  (SST),  are  subgrid-scale  to  the  current  gener¬ 
ation  of  coupled-climate  models,  and  thus  improved 
parameterizations  of  OBLs  with  wave  effects  are  gener¬ 
ally  needed  (Sullivan  and  McWilliams  2010;  Belcher  et  al. 
2012). 

Observations  of  wavy  upper-ocean  dynamics  are  con¬ 
ducted  in  the  presence  of  broadband  wind  wave  com¬ 
plexity,  for  example,  wind  wave  disequilibrium,  persistent 
swell,  high  winds  with  inertially  rotating  wind  stress,  and 
sharp  fronts  and  submesoscale  motions  with  breaking 
waves  are  common  (Smith  1998;  D’Asaro  2001;  Melville 
et  al.  2005;  Capet  et  al.  2008;  Sanford  et  al.  2011).  LES 
studies  of  the  OBL  that  include  vortex  force  are,  how¬ 
ever,  almost  exclusively  conducted  for  idealized  steady 
flow  regimes.  Essentially,  lengthy  LES  time  integrations 
are  carried  out  assuming  a  false  independence  of  winds 
and  waves;  that  is,  there  is  no  wind  wave  evolution  over 
the  long  time  periods  used  to  gather  turbulence  statistics. 
This  is  an  acceptable  approximation  for  near  wind-wave 
equilibrium  conditions  but  is  less  certain  in  nonstationary 
wave  states.  Thus,  a  next  step  in  complexity  is  exploring 
the  dynamics  of  Langmuir  turbulence  in  an  evolving  wind 
wave  regime.  Recent  LES  investigations  consider  wavy 
OBLs  with  a  modest  degree  of  nonequilibrium  wind  wave 
forcing,  for  example,  idealized  wind-wave  misalignment 
(Van  Roekle  et  al.  2012),  nonstationary  mixing  in  a  shal¬ 
low  ocean  (Kukulka  et  al.  2010),  evolution  of  Langmuir 
turbulence  over  a  range  of  wind  speeds  and  wave  age 
(Harcourt  and  D’Asaro  2008),  and  wind-driven  mixing 
below  the  thermocline  (Grant  and  Belcher  2011). 

The  present  study  expands  our  previous  LES  model¬ 
ing  investigations  of  Langmuir  turbulence  to  the  high 
wind  (hurricane)  driven  ocean  boundary  layer.  This  re¬ 
gime  features  rapidly  turning  wind  stress  and  temporally 
evolving  wave  fields  that  are  a  complex  interacting  mix 
of  wind  waves  and  swell  with  multiple  spectral  peaks 
propagating  in  directions  misaligned  with  the  surface 
winds  (Wright  et  al.  2001;  Fan  et  al.  2009).  The  hurricane- 
forced  OBL  then  allows  an  assessment  of  Langmuir 


turbulence  in  a  highly  nonequilibrium  wind  wave  re¬ 
gime.  The  classic  study  of  Price  (1981)  and  more  recently 
Sanford  et  al.  (2011)  emphasize  the  importance  of  in¬ 
ertial  resonance  for  ocean  mixing,  that  is,  wind  stress 
rotating  in  the  same  direction  as  the  underlying  currents. 
These  studies,  however,  do  not  consider  OBL  dynamics 
in  the  presence  of  the  companion  significant  surface  wave 
fields  that  high  winds  generate. 

The  goal  of  the  present  work  is  to  examine  the  re¬ 
sponse  of  the  wavy  OBL  to  realistic  time  varying  wind 
and  wave  forcing  induced  by  an  idealized  moving  hur¬ 
ricane.  The  previous  LES  investigation  most  closely  re¬ 
lated  to  the  present  work  is  that  of  Skyllingstad  et  al. 
(2000)  who  examined  OBLs  driven  by  simplified  mono¬ 
chromatic  wave  fields  and  perfectly  resonant  wind  stress, 
that  is,  a  time- varying  surface  stress  with  period  T  =  lirlf 
in  which /is  the  Coriolis  parameter.  A  preliminary  LES 
of  a  hurricane  OBL  with  a  simplified  time- varying  wave 
field  is  described  by  Sullivan  and  McWilliams  (2010). 
There  are  numerous  studies  of  high-wind  OBLs  carried 
out  using  ID  single-column  parameterizations  (e.g.,  Price 
et  al.  1986;  Large  et  al.  1994;  Crawford  and  Large  1996; 
Zedler  et  al.  2002).  These  investigations  emphasize  the 
importance  of  inertially  rotating  wind  stress  for  mixing 
at  the  thermocline  but  do  not  explicitly  include  surface 
wave  effects;  wave  effects  may  be  indirectly  included  in 
these  models  through  empirically  tuned  constants  and 
mixing  rules.  The  computational  elements  of  the  present 
problem  are  a  spectral  wave  prediction  model  WaveWatch 
III  (Tolman  2002;  Romero  and  Melville  2010b)  and  an 
LES  model  of  the  OBL  with  wave  effects  (Sullivan  et  al. 
2007).  The  manuscript  is  arranged  as  follows:  section  2 
describes  the  large-scale  hurricane  wind  and  wave  fields; 
sections  3  and  4  provide  an  overview  of  the  problem 
posing  and  discuss  the  external  wind  and  wave  forcings 
applied  to  the  LES  model;  sections  5  and  6  are  devoted 
to  a  discussion  of  the  results;  and  section  7  is  a  summary 
of  the  findings. 

2.  Hurricane  wind  and  wave  fields 

Present  computer  power  is  insufficient  to  simulta¬ 
neously  resolve  all  dynamical  scales  in  a  hurricane-driven 
ocean  basin  where  the  largest  scales  of  motion  are  hun¬ 
dreds  of  kilometers  in  horizontal  directions  while  the 
smallest  scales  associated  with  Langmuir  turbulence  are 
about  one  meter.  To  make  our  small-scale  turbulence 
computations  tractable,  we  utilize  a  simplified  one-way 
coupling  strategy  where  the  wave  fields  are  first  com¬ 
puted  over  a  large  ocean  basin  using  surface  winds  that 
are  an  idealization  of  a  translating  hurricane  vortex.  In 
the  present  study,  the  wind  fields  come  from  hindcast 
products  and  are  an  idealization  of  Hurricane  Frances. 
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The  winds  are  first  used  to  drive  the  wave  model,  and 
then  the  Stokes  drift,  computed  from  the  wave  spectra, 
and  the  surface  wind  field  are  imposed  as  time  and  space 
varying  external  point  forcings  to  a  LES  model  of  the 
upper  OBL.  This  one-way  coupling  allows  us  to  examine 
how  the  OBL  dynamics,  including  Langmuir  turbulence, 
respond  to  rapidly  varying  large-scale  forcings.  This  cou¬ 
pling  strategy  with  LES  neglects  some  large-scale  process¬ 
es,  namely,  Ekman  pumping  and  mean-field  advection, 
which  are  most  important  for  the  evolution  of  the  spa¬ 
tially  3D  hurricane  wake  at  late  times  (Price  1981;  Sanford 
et  al.  2011). 

Hurricane  Frances  was  a  category  4  hurricane  that 
developed  in  the  Atlantic  basin  in  2004  and  was  one  of 
the  most  heavily  studied  storms  in  the  Coupled  Boundary 
Layer  Air-Sea  Transfer  (CBLAST)  program  (Black  et  al. 
2007).  Novel  profiling  floats  developed  by  Sanford  et  al. 
(2007),  D’Asaro  et  al.  (2007),  and  Sanford  et  al.  (2011) 
document  the  upper-ocean  mixing  induced  by  Frances 
over  several  days  including  the  development  of  the  hur¬ 
ricane  cold  wake.  The  impact  of  the  storm  on  the  large- 
scale  currents  and  scalars  is  modeled  by  Zedler  (2007) 
and  Sanford  et  al.  (2011). 

a.  Winds  and  surface  fluxes 

Surface  winds,  momentum,  and  heat  flux  parameter- 
izations  are  needed  to  drive  the  LES  and  also  the  wave 
model.  We  make  particular  choices  for  our  modeling  but 
appreciate  that  there  is  uncertainty  in  specifying  these 
external  forcings  owing  to  the  variability  in  the  available 
observations.  An  acceptable  representation  of  the  Hur¬ 
ricane  Frances  winds  is  derived  from  the  NOAA  hindcast 
data  products  described  by  Zedler  (2007,  see  appendix). 
The  specific  formulas  for  the  tangential  and  inflow  radial 
winds  (i/*,  i/)  expressed  in  terms  of  the  radius  r  are 


*  =  PaC<t  |U10|U10>  (2) 

where  pa  is  the  air  density  and  Uio  =  (U1Q,  Vi0)  is  the 
surface  wind  at  a  reference  height  za  -  10  m.  The  drag 
coefficient  is  parameterized  with  the  Large  and  Pond 
(1981)  prescription 

f  0.0012,  |U10|  <llms-1 

d  | (0.49  +  0.065|U10|)  X 10-3,  |U10|  < 25.0m s"1. 

(3) 


This  formula  likely  overpredicts  the  drag  coefficient  at  very 
low  winds  |Ui0|  —  5  m  s_1  but  closely  follows  the  obser¬ 
vations  at  wind  speeds  10  ^  |Ui0|  ^  20  m  s”1.  Based  on 
a  growing  body  of  experimental  data  collected  in  the  field 
and  laboratory  (Powell  et  al.  2003;  Donelan  et  al.  2004; 
French  et  al.  2007)  we  further  limit  (saturate)  the  maximum 
Cd  -  1.8  X  10-3  at  wind  speeds  beyond  20  m  s_1.  This 
saturated  drag  coefficient  is  near  the  estimates  ob¬ 
tained  from  field  observations  but  is  lower  than  the 
value  found  by  Donelan  et  al.  (2004),  Cd  ~  2.3  X  10-3, 
and  Bell  et  al.  (2012,  see  Fig.  19)  who  report  a  mean 
Cd  —  2.4  X  10-3  above  50  m  s_1  based  on  a  budget  anal¬ 
ysis  in  the  atmosphere.  Sanford  et  al.  (2011,  see  Fig.  7), 
using  an  integral  momentum  balance  in  the  ocean,  finds 
smaller  values  Cd  —  1.5  X  10“3  in  the  wind  speed  range 
from  30  to  50  m  s“l. 

Surface  latent  Q\at  and  sensible  Qsen  fluxes  are  also 
imposed  at  the  water  surface  in  the  LES.  The  parame¬ 
terization  of  these  fluxes  is 

e«.=p.c*cp|uI0|Ar. 

e,at=PaCeHU10|A<7,  (4) 


i/(r)  =  53.86e“00054r+2.86; 

i/(r)  =  13.27e~°'0035r— 3.33,  (1) 

which  apply  over  the  domain  rmax  <  r  <  20  rmax,  where 
''max  -  40  km  is  the  radius  of  maximum  winds.  As  sug¬ 
gested  by  Zedler  (2007),  these  expressions  are  aug¬ 
mented  with  a  decaying  exponential  function  for  r  < 
''max-  In  (1),  the  winds  are  in  units  of  meters  per  second, 
r  is  in  units  of  kilometer,  and  the  wind  vr  is  positive  in¬ 
ward.  Based  on  the  measurements  of  Sanford  et  al. 
(2007),  we  also  add  a  constant  translation  speed  Vu  — 
5.5  m  s  ~ 1  to  (1),  which  results  in  an  asymmetrical  storm. 

The  parameterization  of  the  surface  stress  t  at  high 
winds  is  a  key  atmosphere-ocean  coupling  variable  in 
hurricane  prediction.  There  is  vigorous  debate  concern¬ 
ing  its  parameterization  but,  as  is  customary,  we  adopt  the 
bulk  aerodynamic  formula 


where  Cp  is  the  specific  heat,  Lv  is  the  latent  heat  of 
vaporization,  (A  T,  A q)  are  the  temperature  and  specific 
humidity  differences  between  the  sea  surface  and  the 
reference  height,  and  (C*,  Ce )  are  exchange  coefficients. 
We  choose  Ch  -  Ce  -  1.5  X  10-3,  A T  -  2.5  K,  and 
compute  A q  based  on  a  relative  humidity  of  80%  at  the 
reference  height  and  assuming  saturation  at  the  sea  sur¬ 
face.  These  choices  are  representative  values  used  in  hur¬ 
ricane  models  (e.g.,  Price  1981;  Jacob  et  al.  2000;  Jacob 
and  Shay  2003;  D’Asaro  2003;  Zedler  et  al.  2009)  and 
result  in  a  total  heat  flux  from  the  ocean  to  the  atmosphere 
of  more  than  1000  W  m-2  at  the  time  of  maximum  winds. 
Despite  this  seemingly  high  value  the  surface  heat  flux  is 
—0.1  of  the  dominant  entrainment  flux  at  the  time  of 
maximum  winds.  The  sensitivity  of  the  ocean  boundary 
layer  evolution  to  the  surface  exchange  coefficients  is 
described  by  Zedler  et  al.  (2009). 
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Fig.  1.  The  hurricane  storm  propagates  from  left  to  right,  and  the 
fields  are  shown  in  the  form  of  an  X-t  Hovmoller  diagram;  that  is, 
the  temporal  variation  of  a  particular  field  is  shown  for  a  chosen  X 
cross  section.  Time  is  referenced  to  the  time  of  maximum  winds  tm. 
Wind  speed  |Ui0|  (color  scale  m  s  *)  imposed  in  the  WaveWatch 
III  simulations  and  (bottom)  the  wave  age  A  -  CpI |UI0|  over  the 
lifetime  of  the  storm. 


b.  Wave  field 

The  surface  wind  and  wave  fields  input  to  the  LES 
are  products  derived  from  the  wave  prediction  model 
WaveWatch  III  (Tolman  2002)  adapted  to  hurricane 
conditions  (Romero  and  Melville  2010b;  L.  Romero  et  al. 
2012,  manuscript  in  preparation).  The  hurricane  wave  sim¬ 
ulations  use  nearly  exact  computations  of  the  nonlinear 
energy  transfers  (van  Vledder  2006),  wind  input  by  Snyder 
et  al.  (1981),  and  the  energy  dissipation  by  Romero  and 
Melville  (2010b)  with  the  following  modifications.  For 
this  study,  the  dissipation  is  improved  to  allow  a  smoother 
development  of  the  spectrum  with  multiple  peaks  and  to 
allow  the  spectral  tail  to  follow  the  rapidly  turning  winds, 
while  matching  the  empirical  degree  of  saturation  with 
weak  dependence  on  the  wind  speed  or  friction  velocity 
(Banner  et  al.  1989;  Romero  and  Melville  2010a).  A  de¬ 
scription  of  the  code  and  its  validation  for  high  winds  and 
hurricanes  is  given  by  Romero  and  Melville  (2010b)  and 
Romero  et  al.  (2012). 

For  the  present  application,  the  wind-wave  simula¬ 
tions  of  Hurricane  Frances  are  carried  out  in  an  X-Y 
domain  (±410,  ±810)  km  with  a  horizontal  resolution  of 


storm  propagation 
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FIG.  2.  Stokes  velocity  components  (top)  us  and  (bottom)  vr  at  the 
water  surface  for  the  imposed  wind  field  (color  scale  units  ms-1). 


10  km,  that  is,  83  X  163  grid  points.  The  spectral  grid  is 
discretized  with  a  directional  resolution  of  10°  using 
a  constant  relative  bandwidth  dmlm  —  0.078;  w  is  the 
radial  frequency.  The  range  of  resolved  frequencies  is 
[0.157, 15.341]  rad  s-1  corresponding  to  the  wavenumber 
range  [0.0025,  24]  rad  m  \  The  storm  eye  is  initially 
placed  at  (X,  Y)  =  (0,  -510)  km  and  is  translated  in  time. 
Time  series  of  winds  and  waves  are  collected  at  61  dis¬ 
crete  points  distributed  along  the  line  ( X ,  Y)  =  (±300, 0) 
km.  The  primary  output  of  the  wave  model  is  the  2D 
wave  height  spectrum  F(m,  <£),  which  varies  with  w  and 
angular  position  <f> .  The  wave  spectra  are  postprocessed 
to  generate  the  Stokes  drift  velocity  and  wave  age  fields. 

The  X-t  Hovmoller  contour  maps  presented  in  Figs. 
1  and  2  illustrate  some  of  the  variability  in  the  wind  and 
wave  fields  induced  by  our  hypothetical  storm  over  the 
horizontal  domain  X  =  ±300  km  for  a  time  period  of 
more  than  50  h.  In  these  and  subsequent  figures,  time  is 
referenced  to  the  time  of  maximum  winds  tm  on  the 
right-hand  side  of  the  storm.  Here  tm  =  45.625  h  is  based 
on  a  constant  storm  translation  speed  V\x  -  5.5  m  s"1 
and  the  initial  storm  center  location.  Thus,  the  cross  sec¬ 
tion  ( X ,  y  “  0)  is  ahead  of  the  storm  for  t  -  tm<  0,  cuts 
through  the  storm  eye  at  t  -  tm  -  0,  and  is  behind  the 
storm  (in  the  storm  wake)  for  t  -  tm  >  0.  Inspection  of 
the  contours  in  Fig.  1  shows  that  the  wind  speed  [Ujo|  is 
asymmetrical  about  X  =  0  km,  a  consequence  of  the 
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translation  speed  Vh,  with  the  largest  values  —50  m  s”1 
on  the  right-hand  side  of  the  storm  track. 

Wave  age,  shown  in  the  lower  panel  of  Fig.  1,  is  de¬ 
fined  as  A  =  Cp/|Uio]  where  Cp  is  the  phase  speed  of  the 
peak  in  the  wave  height  spectrum:  A  varies  widely  in 
space  and  time  across  the  storm  path.  Extensive  periods 
of  growing  waves  A  —  0.5  are  observed  near  the  time  of 
maximum  winds  with  the  lowest  wave  age  on  the  right- 
hand  side  of  the  storm.  For  time  periods  well  ahead  of 
the  storm  center,  A  —12  indicative  of  near  wind-wave 
equilibrium  conditions.  At  -22  h,  the  wave  age 

jumps  abruptly,  A  >  2,  as  fast-moving  swell  propagates 
rapidly  ahead  of  the  storm.  In  the  wake  of  the  storm  t  - 
tm  >  0,  the  wave  age  again  slowly  approaches  wind  wave 
equilibrium  and  at  long  times  t-  tm>  20  h  the  wave 
state  again  becomes  swell  dominated.  These  estimates  of 
wave  age  in  a  hurricane  are  consistent  with  values  deduced 
from  observations  of  ID  spectra  documented  by  Young 
(1998, 2006)  who  reports  values  ranging  from  [0.5  -  2.0]. 

A  key  variable  for  our  LES  is  the  Stokes  drift  velocity  us , 
which  appears  in  the  vortex  force  u5  X  Z5  where  a>  is  the 
resolved  vorticity;  see  (8).  For  a  broadband  wave  spec¬ 
trum,  ns  is  obtained  by  the  integration  (Kenyon  1969): 

u s(z)  -  2 1  |  kmF(ro,  <f>)e2^z  d<f>  dm,  (5) 


FIG.  3.  (top)  The  turbulent  Langmuir  number  La,  (nondimen- 
sional  color  scale)  and  (bottom)  Stokes  depth  scale  Ds .  In  the 
bottom  panel  the  color  bar  is  in  units  of  meters. 


where  k  denotes  the  horizontal  wavenumber  vector 
[also  see  McWilliams  and  Restrepo  (1999)].  The  Stokes 
drift  velocity  us  inherits  the  complexity  of  the  wave  field 
through  (5);  it  varies  with  wind  speed,  spatial  position, 
and  time  as  depicted  in  Fig.  2.  For  our  particular  choice 
of  storm  propagation  direction,  the  largest  Stokes  drift 
occurs  in  component  vs  on  the  right-hand  side  of  the  storm 
track.  The  distribution  of  the  wave  field  is  more  asym¬ 
metrical  than  the  wind  field  mainly  due  to  the  “extended” 
fetch  on  the  right-hand  side  of  a  translating  hurricane 
(Young  2003).  The  asymmetrical  distribution  of  vs  about 
X  =  0  km  reflects  the  stronger  winds  and  better  wind 
wave  alignment  on  the  storm  right  hand  side.  Clearly  us 
is  largest  (most  negative)  in  front  of  the  storm  center 
with  contours  roughly  symmetrical  about  X  =  0  km. 

McWilliams  et  al.  (1997)  introduced  a  surface  based 
turbulent  Langmuir  number 


La, 


|u,(*.z  =  0,0l 


(6) 


as  a  measure  of  the  relative  importance  of  shear  forc¬ 
ing  to  vortex  force.  In  (6)  the  water  friction  velocity 
u*(Y,  t)  —  y/\r\ !pw,  where  pw  is  the  density  of  water,  u5 
is  evaluated  at  the  water  surface,  and  t  is  given  by  (2). 
Grant  and  Belcher  (2009)  also  adopt  (6)  in  their  scaling 


of  the  turbulent  kinetic  energy  equation  with  wave  ef¬ 
fects.  Harcourt  and  D’Asaro  (2008)  prefer  a  surface  based 
Langmuir  number  LasL  =  yju*i{u 5)SL,  where  {w*)s l  is  the 
average  Stokes  drift  in  the  surface  layer.  These  definitions 
of  Langmuir  number  are,  however,  strictly  only  applicable 
in  the  situation  of  aligned  winds  and  waves;  Van  Roekle 
et  al.  (2012)  attempts  to  generalize  (6)  for  wind-wave 
misalignment.  Here  we  adopt  (6)  as  our  working  defini¬ 
tion  of  Langmuir  number  because  of  its  simplicity  and  the 
ability  to  compare  with  previous  studies;  note  that  in  (6)  us 
is  computed  from  the  full  wave  spectrum  according  to  (5). 

Figure  3  shows  contours  of  turbulent  Langmuir 
number  for  our  hypothetical  storm.  Outside  the  radius 
of  maximum  winds  there  are  broad  regions  where  0.25  < 
La,  <  0.4.  The  very  low  values  in  the  storm  eye  should  be 
favorable  for  generating  Langmuir  cells  with  high  aspect 
ratio,  i.e.,  with  their  streamwise-to-crosswise  length  scale 
much  greater  than  unity  (Chini  et  al.  2009).  Historically 
the  criterion  La,  <  0.4  (e.g.,  Li  et  al.  2005),  used  for 
aligned  winds  and  waves,  implies  that  often  the  hurri¬ 
cane  driven  OBL  should  be  in  a  regime  favorable  to  the 
generation  of  Langmuir  turbulence.  However,  the  ap¬ 
plicability  of  the  criterion  La,  <  0.4  for  nonequilibrium 
situations  where  weaker  Langmuir  cells  can  develop  off- 
axis  between  the  wind  and  wave  directions  is  unknown 
(Van  Roekle  et  al.  2012). 
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The  vertical  profile  of  Stokes  drift  uy(z)  spirals  with 
depth,  a  consequence  of  a  multicomponent  wave  field 
composed  of  swell  and  wind  waves  that  propagate  in 
different  directions.  To  characterize  the  depth  penetra¬ 
tion  by  the  Stokes  drift,  an  important  variable  in  Lang¬ 
muir  turbulence,  we  define  a  vertical  length  scale 


■  J,  iu 


0  •«,(<>.  0 


«>,or 


dz. 


(7) 


where  the  integral  is  taken  over  the  interval  zs  <  z  <  0 
with  zs  chosen  as  the  depth  where  the  integrand  remains 
positive  [see  Harcourt  and  D’Asaro  (2008)  for  an  al¬ 
ternate  length  scale  for  aligned  winds  and  waves].  The 
above  definition  is  robust  and  reduces  to  Ds  —  A  =  2ir!k 
for  a  monochromatic  wave  field.  The  variation  of  DS{X ,  t) 
is  given  in  Fig.  3.  Over  the  time  interval  (-20, 10)  h,  Ds  ~ 
60  m  with  values  greater  than  130  m  on  the  right-hand 
side  of  the  storm  at  the  time  of  maximum  winds.  At  t  —  tm 
5  h,  there  is  a  small  region  to  the  left  of  the  storm  eye 
where  Ds  becomes  small,  a  reflection  of  the  confused 
sea  state  induced  by  opposing  winds  and  waves.  This  is 
also  a  region  where  La,  >  0.6.  Based  on  the  results  in  Fig.  3 
we  conclude  that  the  vertical  profile  of  Stokes  drift  is 
track  dependent  and  differs  noticeably  from  a  monochro¬ 
matic  prescription,  but  occupies  a  significant  fraction  of 
the  OBL  even  at  the  time  of  maximum  winds;  that  is, 
DJ\h\  >  1  where  h  is  the  OBL  depth. 


point,  second-moment  turbulent  kinetic  energy  (TKE) 
closure  subgrid-scale  parameterization  and  a  flat  ocean 
surface  (Moeng  1984;  McWilliams  et  al.  1997).  The  added 
wave  effects  are  the  vortex  force,  Stokes-Coriolis  cou¬ 
pling,  Lagrangian  mean  advection  associated  with  Stokes 
drift,  a  wave-averaged  increment  to  the  pressure  that 
arises  through  conservative  wave-current  interaction 
(McWilliams  et  al.  1997),  and  Stokes  production  in  the 
subgrid-scale  (SGS)  TKE  equation  (Sullivan  et  al.  2007). 
Additional  acceleration  and  energy  generation  due  to 
nonconservative  stochastic  wave  breaking,  described  by 
Sullivan  et  al.  (2004,  2007),  are  not  considered  here. 

We  use  standard  LES  notation  where  spatially  filtered 
(or  resolved)  variables  are  denoted  with  an  overbar: 
(u,w  =  Vxn)  are  the  resolved-scale  velocity  and  vortic- 
ity;  e  is  the  subgrid-scale  TKE;  0  is  the  temperature;  and 
7r  is  a  generalized  pressure  variable.  Other  variables  ap¬ 
pearing  in  the  LES  equation  set  are  the  Coriolis  frequency 
f  =  (0, 0,/);  reference  temperature  Q0\  gravitational  ac¬ 
celeration  g;  and  subgrid-scale  momentum  and  density 
fluxes,  (r,y,  r,e),  respectively.  These  subgrid-scale  fluxes 
are  modeled  using  the  eddy  viscosity  prescription  de¬ 
scribed  in  Moeng  (1984)  and  Sullivan  et  al.  (1994),  which 
implies  that  the  principal  feedback  of  e  on  u  and  0  occurs 
through  subgrid-scale  mixing  with  eddy  viscosity  and 
diffusivities  (ph  vs)  oc  el/2. 

In  the  situation  of  misaligned  winds  and  waves,  Stokes 
drift  velocity  couples  to  the  resolved  vorticity  =  (£,  17,  £) 
to  generate  a  three-dimensional  vortex  force 


3.  LES  problem  posing 

The  structure  of  the  hurricane  driven  OBL  depends 
on  the  large-scale  forcings  discussed  in  section  2.  To 
thoroughly  investigate  all  possible  combinations  of  wind 
and  wave  forcing  is  beyond  our  computational  reach 
with  the  LES  and,  thus,  we  restrict  our  attention  to  two 
locations  X  —  ±60  km  from  the  storm  center  shown  in 
Fig.  1.  These  locations  are  of  interest  because  they  fea¬ 
ture  high  wind  conditions  with  different  wind  stress  time 
histories  and  have  a  mixture  of  differently  oriented  winds 
and  waves.  As  a  further  parameter  variation,  on  each  side 
of  the  storm,  simulations  are  also  conducted  with  and 
without  surface  wave  effects,  that  is,  with  and  without 
Stokes  drift  which  generates  the  vortex  force. 

a .  LES  model 

The  details  of  our  LES  model  for  the  OBL  with  wave 
effects  are  extensively  discussed  by  McWilliams  et  al. 
(1997)  and  Sullivan  et  al.  (2007).  The  equations  for  the 
resolved  flow  components  are  the  Craik-Leibovich  (CL) 
theory  with  crucial  wave-current  coupling  through  the 
vortex  force.  The  LES  model,  with  wave  effects,  is  based 
on  the  incompressible  Boussinesq  equations  with  a  single¬ 


u,x“  =  (vsC,  -usC,  us7]  -  Vs£)  (8) 

that  contributes  to  all  three  resolved  momentum  equa¬ 
tions;  in  the  case  of  aligned  winds  and  waves,  which  is 
most  often  studied,  the  vortex  force  reduces  to  the  two- 
component  form  1^(0,  —  f,  17). 

b .  LES  setup 

Two  computational  domains  of  varying  size  are  used 
in  the  LES.  First,  the  vertical  dimension  of  the  compu¬ 
tational  box  is  chosen  to  be  about  twice  the  expected 
maximum  depth  of  the  OBL;  based  on  observations 
l^lmax  ~  120  m  (Sanford  et  al.  2011)  and  thus  we  choose 
a  box  of  depth  —240  m.  The  horizontal  dimension  is 
picked  to  be  — 5|/i|max  so  as  to  allow  the  interior  motions 
to  develop  independent  of  the  periodic  sidewall  boundary 
conditions  (see  Schmidt  and  Schumann  1989).  On  the  left- 
hand  side  X  =  —  60  km,  which  has  slightly  weaker  forcing, 
a  small  domain  ( Lx ,  Ly,  Lz)  —  (750,  750,  —240)  m  was 
deemed  acceptable.  On  the  right-hand  side  of  the  storm, 
X  -  60  km,  the  vigorous  forcing  from  the  inertially  ro¬ 
tating  wind  stress  was  found  to  generate  a  system  of 
large-scale  random  gravity  waves  in  the  hurricane  wake. 
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Test  runs  indicated  that  these  structures  were  only  ad¬ 
equately  captured  by  increasing  the  horizontal  area  of 
the  small  domain  LES  by  a  factor  of  4.  Thus,  a  large 
domain  is  used  on  the  right  hand  side  with  ( Lx ,  Ly ,  Lz)  = 
(1500, 1500,  -240)  m.  The  number  of  grid  points  in  the 
small  domain  is  ( Nx ,  Nyt  Nz)  =  (512,  512,  256)  and  is 
increased  to  (1024, 1024,  256)  in  the  large  domain;  thus 
the  horizontal  spacing  Ajc  =  A y  -  1.46  m  in  all  the  com¬ 
putational  boxes.  In  the  vertical  direction,  we  use  a 
stretched  vertical  grid  where  the  spacing  between  any 
neighboring  cells  is  held  fixed  at  the  ratio  K  -  1.00451 
with  the  first  w  grid  point  located  at  Azi  =  -0.5  m;  thus 
the  first  u  point  is  at  z\i 2  =  -0.25  m  in  a  staggered  grid. 
Based  on  these  choices,  the  cell  aspect  ratio  AxIAz  = 
(2.93,  1.54)  at  z  -  (0,  -100)  m.  The  small  and  large 
computational  boxes  are  able  to  capture  small-scale 
turbulence  at  early  times  and  larger  scales  of  motion  that 
develop  during  the  rapid  deepening  of  the  OBL  at  the 
time  of  maximum  winds. 

The  initial  temperature  sounding  is  simply  con¬ 
structed  from  two  linear  segments:  over  the  depth  hi  ^ 
z  <  0,  Qi  =  302.4  K  and  for  z  <  ht  the  temperature  de¬ 
creases  at  a  rate  of  0.04  K  m-1.  This  is  approximately 
the  variation  observed  for  Hurricane  Frances  (Sanford 
et  al.  2011).  The  depth  ht  -  -32  m  is  the  initial  OBL 
depth.  The  Coriolis  parameter  is  /  -  0.681  X  10-4  s-1, 
corresponding  to  an  inertial  period  T  —  25.6  h.  The 
simulations  are  initiated  from  rest  with  small  random 
temperature  perturbations  under  the  wind  and  wave 
forcing  from  hour  10  of  the  WaveWatch  III  simulation. 
At  this  time,  the  wind  and  wave  fields  are  slowly  varying 
with  winds  less  than  10  m  s"1.  Next,  these  external  forc¬ 
ings  are  held  constant  and  the  simulations  are  inte¬ 
grated  for  10  physical  hours  to  generate  fully  developed 
equilibrium  turbulence.  The  resulting  currents  are 
—  0.05  ms1,  and  the  initial  state  of  these  weak  currents 
has  minimal  impact  on  the  long  time  ocean  response 
under  strong  forcing  (Crawford  and  Large  1996,  p.  889). 
Finally,  the  winds  and  Stokes  drift  are  allowed  to  follow 
the  time  history  shown  in  Figs.  5  and  6.  The  time  step  At 
in  the  code  is  dynamically  picked  to  obey  a  fixed  Courant- 
Fredrichs-Lewy  condition.  Here  At  varies  by  an  order  of 
magnitude,  from  0.3  to  5  s,  and  more  than  300  000  time 
steps  are  needed  to  cover  the  70  physical  hours  of  the 
simulation.  With  a  storm  translation  speed  of  5.5  m  s”1 
the  LES  domains  are  initially  located  more  than  900  km 
from  the  storm  center.  The  parallelization  of  the  code  is 
described  by  Sullivan  and  Patton  (2011). 

4.  Transient  wind  and  wave  forcing  applied  to  LES 

The  temporal  and  spatial  variability  in  the  large-scale 
forcing  conditions  critically  impacts  the  time  evolution 


FIG.  4.  External  winds  (Uio,  Pio)  and  cooling  Q*  applied  to  the 
LES:  (a)  right-hand  side  of  the  storm  X  =  60  km  and  (b)  left-hand 
side  of  the  storm  X  =  -60  km. 


of  the  OBL  and  the  development  of  Langmuir  turbu¬ 
lence.  To  aid  the  interpretation  of  the  OBL  dynamics 
generated  by  the  simulations  we  first  discuss  the  time 
history  of  winds,  waves,  and  cooling  that  are  applied  at 
the  upper  boundary  of  the  LES;  see  Figs.  4, 5, 6,  and  7.  In 
these  figures,  the  winds  Uio,  surface  cooling  flux  Q*  = 
2sen  +  2iat>  Stokes  drift  velocity  u5,  wind-wave  align¬ 
ment,  and  surface  wind  stress  t  are  time  slices  at  the  hor¬ 
izontal  locations  X  =  ±60  km  taken  from  Figs.  1  and  2. 

The  primary  forcing  to  the  LES  (and  also  the  wave 
model)  are  the  surface  winds,  shown  in  Fig.  4,  converted 
to  a  surface  wind  stress  according  to  the  prescription  (2). 
The  amplitude  of  the  wind  stress,  depicted  in  Figs.  5  and 
6,  is  a  strong  function  of  (. X ,  t)  within  the  storm  track  and, 
as  expected,  each  component  displays  a  strong  pulse  at 
the  time  of  maximum  winds  mimicking  the  surface  winds. 
Figures  5  and  6  also  compare  the  time  varying  stress  t(*) 
with  an  “inertial”  stress  of  equal  magnitude  that  rotates 
clockwise  with  frequency/,  that  is,  t,(/)  =  |T(/)|e";^.  This 
exposes  the  direction  and  rate  of  rotation  relative  to  an 
inertial  motion.  Both  wind  stress  components  on  the  right 
side  of  the  storm  track  are  partially  resonant  with  the 
ideal  inertial  stress  t,.  The  period  of  resonance  is  rela¬ 
tively  short,  about  10  h,  but  occurs  during  the  time  of 
maximum  winds  (-5,  5)  h.  Meanwhile,  on  the  left-hand 
side,  the  wind  stress  magnitudes  are  modestly  lower  than 
on  the  right-hand  side  and  component  ry  clearly  rotates  in 
a  direction  opposite  to  that  of  an  inertial  rotation.  Tran¬ 
sient  wind  forcing  at  frequency  /  with  clockwise  rotation 
is  able  to  drive  currents  below  the  mixed  layer  and  thereby 
efficiently  mix  the  OBL  (e.g.,  Price  1981;  Crawford  and 
Large  1996;  Skyllingstad  et  al.  1999).  The  observations 
of  Sanford  et  al.  (2011)  demonstrate  how  wind  stress 
with  limited  resonance  time  is  able  to  efficiently  mix 
a  hurricane  OBL.  At  early  and,  more  so,  at  late  times  the 
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FlG.  5.  Surface  wind  stress  and  wave  forcing  on  the  right  hand 
side  of  the  storm,  X  =  60  km:  (a)  time  variation  of  the  wind  stress  t 
where  the  (black,  gray)  lines  are  the  components  r  -  (t*,  r^),  re¬ 
spectively.  The  (black,  gray)  dotted  lines  are  resonant  stress  com¬ 
ponents  T/(t)  that  rotate  clockwise  with  inertial  frequency/;  that  is, 
T|(f)  =  \T(t)\e~ift.  (b)  The  surface  Stokes  drift  velocity  components 
(f/j,  w,)  normalized  by  the  local  friction  velocity  u+(X,  t)  are  de¬ 
picted.  (c)  The  wind-wave  alignment  given  by  (9)  at  the  water 
surface  (black  line)  and  at  depth  z  =  -(5.8, 12.4, 40,  80.2)  m  gray 
(solid,  dotted,  dash-dot,  long  dash)  lines,  respectively. 


wind  stress  rotation  applied  at  the  upper  boundary  of  the 
LES  is  noninertial  on  both  sides  of  the  storm  track. 

Surface  cooling  plays  a  minor  role,  despite  the  more 
than  1200  W  m~2  applied  at  the  time  of  maximum  winds; 
it  is  offset  by  the  rapid  increase  in  surface  shear  stress. 
The  ratio  \hj\IL  where  L  -  —u*IPkQ*  is  the  Monin- 
Obukhov  length  is  often  used  in  atmospheric  boundary 
layer  analysis  as  a  measure  of  the  relative  importance  of 
shear  and  buoyancy  forcing  at  the  surface  (e.g.,  Deardorff 
1972;  Moeng  and  Sullivan  1994).  In  this  expression,  the 
buoyancy  parameter  p  —  glda  and  k  =  0.4  is  the  von 
Karman  constant.  Broadly,  over  the  bulk  of  the  time 
period  the  OBL  surface  shear  forcing  dominates  surface 
cooling  since  -1  <  |/i,*|/L  <  0  and  during  the  time  period 
of  maximum  winds,  (—10, 10)  h,  we  find  —0.1  <  \hj\IL  <  0. 
Cooling  plays  a  modest  role  at  early  times  and  con¬ 
tributes  to  the  well-mixed  (vertically  uniform)  shape  of 
the  mean  current  profiles. 

It  is  informative  to  discuss  how  the  Stokes  velocity 
varies  in  our  simulations  since  it  impacts  the  development 
of  Langmuir  turbulence.  In  the  lower  and  middle  panels 
of  Figs.  5  and  6  we  show  the  time  variation  of  the  wind- 
Stokes  alignment 


cos  0  = 


uio(0-us(z’0 

|U10(OI|u5(z,OI 


(9) 


Fig.  6.  Surface  wind  stress  and  wave  forcing  on  the  left-hand  side 
of  the  storm,  X  =  -60  km.  The  panels,  labeling  and  variables  are 
the  same  as  in  Fig.  5.  Note  the  difference  in  the  vertical  scale  in  the 
lower  panels  of  Figs.  6  and  5. 

at  depths  z  =  (0,  -5.8,  -12.4,  -40,  -80.7)  m,  and  the 
components  of  the  normalized  Stokes  drift  velocity  at 
the  surface  ( us ,  vs)!u+.  The  surface  based  La„  given  by 
(6),  can  be  inferred  from  this  plot  and  its  range  is  narrow; 
along  X  =  60  km,  0.28  <  La,  <  0.32  and  along  X  = 
-60  km,  0.28  <  La,  <  0.41.  First,  we  notice  ny(0,  t)lu* 
has  the  same  shape  and  amplitude  on  both  sides  of  the 
storm  track,  X  =  ±60  km.  This  is  expected  based  on 
the  contours  of  Stokes  drift  velocity  in  Fig.  2.  Thus,  the 
important  differences  in  uju *  between  the  right-hand 
and  left-hand  sides  are  due  to  component  vju *.  On  the 
right-hand  side  of  the  track  Uy(0,  f)/w*  is  nearly  constant 
with  time,  less  than  a  25%  variation  over  the  period 
(—30,  5)  h.  This  correlation  between  Stokes  drift  and 
wind  stress  is  meaningful  since  u *  varies  by  a  factor  of  7 
over  the  time  period  of  the  simulations.  Meanwhile  on 
the  left-hand  side  u5(0,  t)lu*  at  first  decreases  linearly, 
then  changes  sign  and  finally  reaches  a  minimum  near 
t  —  tm  =  0  h.  These  differences  are  a  result  of  the  ex¬ 
tended  fetch  on  the  right  side  of  the  storm  versus  the 
progressive  misalignment  between  winds  and  waves  on 
the  left  side  as  time  advances. 

The  wind-Stokes  alignment  is  higher  at  the  surface 
than  at  depth,  which  results  from  a  multimode  surface 
wave  spectrum  composed  of  wind  waves  and  swell:  re¬ 
call  in  wind  wave  equilibrium  cos <f>  =  1  at  all  z.  Short 
wavelength  wind  waves  rapidly  align  themselves  to 
varying  surface  shear  forcing  compared  to  the  slow  ad¬ 
justment  of  long  wavelength  swell;  this  pattern  of  wave 
development  in  a  hurricane  is  reported  by  Wright  et  al. 
(2001).  Thus  surface  values  of  Stokes  drift  are  only  a 
partial  reflection  of  the  wave  state  that  drives  Langmuir 
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FIG.  7.  Vertical  profiles  of  Stokes  drift  velocity  normalized  by  the  local  friction  velocity  (m„  vs)/u+  at  time  instances 
t  -  tm  -  (-25,  -10,0,5)  h:  z  is  made  dimensionless  by  the  Stokes  depth  scale  Ds  given  by  (7).  The  solid  (black,  gray) 
lines  are  (uSy  vs)iu *,  respectively,  along  the  line  X  -  60  km.  The  dashed  (black,  gray)  lines  are  their  counterparts 
along  the  line  X  =  -60  km.  For  the  selected  time  instances  the  depth  scale  Ds  =  (14.5, 86.0, 123.5, 74.7)  m  and  (12.6, 
75.2, 45.7, 37.1)  m  for  X  -  60  (-60)  km.  Note  at  t  -  tm  =  5  h  the  solid  black  and  gray  lines  overlap. 


turbulence.  Swell  first  appears  at  t  —  tm  —  —22  h,  as  in¬ 
dicated  by  wave  age  A  >  1.2  in  the  lower  panel  of  Fig.  1, 
and  reduces  the  wind-Stokes  alignment  at  depth  through¬ 
out  the  passage  of  the  storm.  At  the  surface,  on  the  right- 
hand  side  of  the  storm  the  wind-Stokes  alignment  is 
always  within  10°  at  the  surface  and  no  more  than  45°  at 
z  —  -5.8  m.  At  the  time  of  maximum  winds  the  Stokes 
drift  is  near  “resonant”  with  the  winds;  that  is,  the  wind- 
Stokes  alignment  is  near  unity  over  a  depth  of  more  than 
80  m  and  the  Stokes  depth  scale  exceeds  130  m.  Thus  we 
anticipate  strong  Langmuir  turbulence  to  develop  at  this 
time.  Over  the  time  period  (-5,  5)  h,  the  wind-Stokes 
alignment  on  the  left-hand  side  increases,  60°  at  the  sur¬ 
face  and  almost  90°  at  z  =  - 10.2  m,  because  the  wind  and 
dominant  wave  propagation  directions  are  opposed. 

Selected  vertical  profiles  of  the  Stokes  drift  velocity 
are  shown  in  Fig.  7  at  times  before,  at,  and  after  the  time 
of  maximum  winds.  These  profiles  illustrate  that  the 
normalized  surface  values  of  uju *  and,  hence,  La„  have 
only  a  mild  variation  over  the  time  period.  The  invariance 
of  La,  however  disguises  the  profile  shape  changes  that 
occur  with  {X,  /).  For  example,  the  profiles  of  vs(z,  t)  at 


t  —  tm-  (—25, 0)  h  have  similar  surface  values  of  La,  but 
the  profile  at  t  -  tm  =  0  extends  much  deeper  into  the 
OBL:  at  t  -  tm  =  (-25,  0,  0)  h  Ds  -  (14.5, 123)  m,  re¬ 
spectively.  Intuitively,  we  expect  the  formation  of  Lang¬ 
muir  turbulence  to  depend  on  La„  the  vertical  shape  of 
the  Stokes  profile  and  its  depth  penetration. 

5.  Identification  of  Langmuir  turbulence 

a.  Temporal  evolution  of  vertical  velocity 

The  LES  of  steady  aligned  wind  wave  flows  and  a  few 
field  observations  find  vertical  velocity  is  a  robust  sig¬ 
nature  of  Langmuir  turbulence.  In  the  Langmuir  tur¬ 
bulence  regime,  coherent  Langmuir  cells  form  and  these 
structures  can  be  detected  in  the  instantaneous  vertical 
velocity  field  and  furthermore  impact  the  statistical  mo¬ 
ments  of  w.  Figure  8  compares  the  maximum  total  (re¬ 
solved  plus  SGS)  second-order  moment  (w2)max  and 
the  minimum  resolved  third-order  moment  (vv3)min  di¬ 
agnosed  from  our  simulations.  Since  the  LES  equations 
are  closed  at  the  second-moment  level  there  is  no 
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FIG.  8.  Temporal  variation  of  the  maximum  total  vertical  velocity  variance  (w2)mux  and 
minimum  resolved  third-order  moment  (vv3)™},,  for  cases  with  and  without  wave  effects  at  lo¬ 
cations  left  and  right  of  the  storm  center.  Simulations  with  wave  effects  for  (a)  X  =  60  km  and 
(b)  X  —  -60  km  and  simulations  with  no  wave  effects  for  (c)  X  =  60  km  and  (d)  X  =  —60  km. 
The  moments  are  normalized  appropriately  by  either  «*(A\f)  or  ul(Xyt). 


estimate  of  third-order  SGS  quantities.  Both  moments 
are  normalized  by  local  u+(X,  /).  These  temporally 
varying  moments  are  obtained  by^  first  computing  the 
vertical  profiles  of  resolved  (w(z,  /)  )  and  { w(z ,  0  )  at  all 
z  for  a  given  /:  angle  brackets  denote  an  x  -  y  spa- 
tial  average  over  the  LES  domain.  The  MO  )max  and 
MO  )min  are  found  by  searching  for  the  appropriate 
maximum  or  minimum  along  their  respective  vertical 
profiles.  We  then  add  an  estimate  of  the  SGS  contribu¬ 
tion,  2e/3,  to  the  resolved  w  variance;  that  is,  <w2)max 
is  the  sum  of  resolved  and  SGS  contributions.  We  note 
that  the  vertical  locations  of  (w2)max  and  <w'3)inin  are  not 
the  same  and  vary  with  time.  These  positions  are  always 
located  in  the  upper  surface  layer  of  the  OBL  but  are 
found  deeper  in  the  water  as  the  winds  and  waves  increase. 

The  time  series  depicted  in  Fig.  8  are  compelling  ev¬ 
idence  for  Langmuir  turbulence  in  the  hurricane  driven 
OBL.  The  w  moments  are  sensitive  to  the  location 
within  the  storm  track,  depend  on  the  time  and  space 
evolving  wind  and  wave  fields,  and  the  inclusion  of  phase 
averaged  wave-current  interactions  in  the  LES.  Over 
the  lifetime  of  the  vortex-force  simulation  (w2)max  >  2.5 
along  the  track  X  =  60  km.  Notice  also  that  during  the 
time  period  with  resonant  wind  forcing,  (—5,  5)  h,  the 
maximum  vertical  velocity  variance  is  near  3  and  this 
level  persists  in  the  hurricane  wake.  Recall  w*  varies  by 


about  a  factor  of  7  and  Ds  varies  from  20  to  130  m  over 
this  time  period  and,  thus,  the  dynamic  range  of  (w2)  is 
broad.  Surprisingly,  simple  scaling  by  w*  captures  much 
of  the  variation  in  w  variance.  On  the  left  side  of  the 
storm  track,  the  maximum  w  variance  at  first  tracks  its 
counterpart  on  the  right-hand  side,  then  gradually  de¬ 
cays,  reaches  a  minimum  value  at  the  time  of  maximum 
winds,  and  gradually  recovers  in  the  hurricane  wake.  In 
the  simulations  without  vortex  force,  the  maximum  w 
variance  is  nearly  independent  of  the  track  position  and 
is  similar  to  what  is  expected  for  a  neutral  flat-wall 
boundary  layer  flow. 

The  differences  between  the  variations  on  the  right 
and  left  hand  side  of  the  storm  can  be  traced  to  the  wave 
fields,  the  wind-Stokes  alignment  shown  in  Figs.  5  and  6, 
and  the  production  of  w  variance  by  Stokes  drift.  The 
Stokes  contribution  to  the  ensemble  average  w  variance 
budget  consists  of  two  terms 

Vs  -  T#*  +  7f)  =  -  (u'w)  ^  -  (v'w)  (10) 

where  {u'w)  and  {u'w)  are  ensemble  average  (total) 
vertical  momentum  fluxes.  Note  there  is  no  mean  shear 
production  in  the  w  budget.  The  derivation  of  (10)  fol¬ 
lows  the  same  sequence  of  steps  used  to  derive  the  TKE 
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FlG.  9.  Production  of  vertical  velocity  variance  by  Stokes  drift. 
The  terms  =  -{ u'w)dusldz  and  =  ~(v'w)dvsldz  are 
computed  at  the  vertical  location  of  maximum  vertical  velocity 
variance  along  the  right-hand  side  (solid  lines)  and  left-hand  side 
(dotted  lines)  of  the  storm  track:  are  indicated  by  black 

(gray)  line  colors.  The  terms  are  made  dimensionless  by  \hj  \ !u\  and 
are  obtained  by  spatial  averaging  over  an  jc  -  y  plane  and  by  time 
averaging  over  45  min.  The  momentum  fluxes,  which  appear  in 
(/*“\  /*u}),  are  the  sum  of  the  LES  resolved  and  SGS  contributions. 
Note  the  periods  where  (7^M\7^u))  are  negative  on  the  left-hand 
side  of  the  storm  track  indicating  that  the  turbulent  fluxes  and 
gradients  of  Stokes  drift  are  countergradient. 

budget  with  wave  effects  (e.g.,  McWilliams  et  al.  1997).  The 
total  TKE  budget  depends  on  the  vertical  gradients  of  the 
currents  and  Stokes  drift,  i.e.,  the  mean  Lagrangian  shear 

cs=^n+us)i+^+vs)y-  (u) 

Notice  that  in  (10)  the  production  term  Vs  depends  on 
the  correlation  between  turbulent  fluxes  and  Stokes  drift 
gradients.  These  terms  can  be  positively  correlated  as  in 
the  case  of  wind  wave  equilibrium,  uncorrelated  or  even 
negatively  correlated:  that  is,  the  momentum  fluxes  can 
be  countergradient  to  the  Stokes  drift  gradient.  A  time 
series  of  Stokes  production,  given  by  (10),  at  the  location 
of  the  maximum  vertical  velocity  variance  is  shown  in 
Fig.  9.  There  is  a  significant  temporal  correlation  be¬ 
tween  the  enhancement  and  reduction  in  w  variance, 
shown  in  Fig.  8,  with  the  time  variation  of  Stokes  pro¬ 
duction.  In  particular  the  variance  on  the  right  and 
left  hand  sides  of  the  storm  track  first  begin  to  depart  at 
/  —  =  — 12  h.  At  this  time  on  the  left-hand  side  of  the 

storm,  the  production  term  reaches  a  positive  maxi¬ 
mum  and  switches  sign  and  becomes  a  destruction 
term.  An  even  steeper  reduction  in  w  variance  occurs  at 
the  time  of  maximum  winds  where  V ^  is  a  sink  of  w 
variance  and  is  small  and  positive.  On  the  right- 
hand  side  of  the  storm  are  always  positive  with 
reaching  a  maximum  in  the  storm  wake.  Thus  Stokes 
production  is  most  often  a  significant  source  of  w  variance, 
but  can  in  certain  misaligned  wind  wave  states  be  a  sink. 


Wave-current  interactions,  which  are  key  processes  in  the 
budgets  of  turbulence,  play  a  role  in  the  enhancement  and 
reduction  of  w  variance  in  the  hurricane  OBL.  The  results 
emphasize  the  importance  of  the  vertical  gradients  of  the 
Stokes  drift  profile  and  their  alignment  with  the  vertical 
momentum  fluxes.  It  is  worth  noting  that  <  0  also 
implies  that  the  currents  and  waves  are  moving  in  op¬ 
posite  directions  near  the  surface  (see  section  6b),  which 
can  enhance  wave  breaking  (Melville  1996). 

D’Asaro  (2001),  using  drifting  Lagrangian  floats,  re¬ 
ports  measurements  of  vertical  velocity  variance  that  are 
“remarkably  well  parameterized  by  wind  stress”  with 
values  ranging  from  1.75  to  2  times  that  measured  in 
solid-wall  turbulent  layers  driven  by  shear.  Our  vortex- 
force  simulations  on  the  right  hand  side  of  the  storm 
track  are  in  good  agreement  with  these  measurements, 
and  we  speculate  that  the  observed  high  correlation  with 
wind  stress  is  a  consequence  of  the  tight  relationship 
between  surface  Stokes  drift  and  w*  for  approximately 
aligned  winds  and  waves,  as  shown  in  Fig.  5.  Harcourt 
and  D’Asaro  (2008)  and  Grant  and  Belcher  (2009)  propose 
a  velocity  scale  for  Langmuir  turbulence,  wL  =  u*/Laf3. 
However  our  range  of  La,  is  somewhat  limited  along 
X  =  60  km  because  of  the  tight  correlation  between  uT 
and  w*;  hence  along  this  track  wL  ~  w*.  We  mention  that 
the  applied  wind  forcing  in  our  simulations  varies  smoothly 
with  time  (there  is  no  assumption  of  wind  wave  equi¬ 
librium)  and  then  the  wave  fields  adjust  naturally  to  the 
changing  local  wind  conditions.  Hence  the  Stokes  drift 
velocity  imposed  in  the  LES  also  varies  smoothly  with 
time  and  changing  winds. 

The  impact  of  Langmuir  turbulence  on  the  third-order 
vertical  velocity  moment  is  impressive.  Figure  8  depicts 
the  minimum  (w3)min  based  on  the  resolved  vertical 
velocity.  With  vortex  force  the  minimum  third-order 
moment  can  be  10  times  larger  than  its  no  vortex  force 
counterpart  at  the  time  of  maximum  winds.  Again  we 
notice  a  significant  reduction  in  (w3)min  on  the  left-hand 
side  of  the  storm  in  the  simulation  that  includes  the  vortex 
force  mirroring  the  reduction  in  (w>2)max.  The  third-order 
moment  of  vertical  velocity  indicates  an  abundance  of 
Langmuir  cells,  and  is  important  as  it  appears  in  the 
turbulent  kinetic  energy  equation.  The  changes  in  (Tv3) 
indicate  that  Langmuir  turbulence  has  modified  the 
vertical  transport,  the  divergence  d{Ww2)ldz ,  in  the  OBL 
in  significant  ways  (Grant  and  Belcher  2009).  We  find 
that  the  moments  (w2)max  and  {w3)rnjn  are  stronger  and 
clearer  indicators  of  Langmuir  circulations  than  the 
nondimensional  skewness  Sw  =  (h>3)/(h>2)3/2,  which  is  of¬ 
ten  used  to  describe  Langmuir  turbulence  (e.g.,  Li  and 
Garrett  1993;  McWilliams  et  al.  1997;  Sullivan  et  al.  2007). 
In  our  simulation  with  vortex  force  on  the  right  hand  side 
of  the  storm  minimum  skewness  varies  over  the  range 
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-1.1  <  Sw  <  —0.6,  which  indicates  that  the  downwelling 
velocities  are  narrower  and  stronger  compared  to  the 
broader  upward  motions. 

6.  Flow  structures 

Flow  visualization  of  the  vertical  velocity  field  pro¬ 
vides  insight  into  the  formation  and  evolution  of  the 
coherent  structures  in  a  high-wind  hurricane  driven  OBL. 
Extensive  interrogation  of  the  LES  volumes  shows  that 
elongated  streaky  structures,  that  is,  coherent  Langmuir 
cells,  are  pervasive  at  all  stages  of  the  storm  evolution. 
The  elliptical  shaped  cells  meander  randomly  in  their 
long-axis  direction  with  occasional  merges  at  downstream 
junctions.  Figure  10  is  a  typical  sample  that  illustrates 
these  processes.  This  image  compares  the  instantaneous 
normalized  vertical  velocity  vv/m*  at  the  time  of  maxi¬ 
mum  winds  on  the  left  and  right  hand  sides  of  the  storm. 
Also,  included  in  this  image  are  the  directions  of  the 
surface  and  Stokes  drift  velocity,  which  provide  a  refer¬ 
ence  for  the  cell  alignment.  To  allow  a  fair  comparison 
of  the  spatial  scales,  the  visualization  is  depicted  over 
equal  subareas  of  the  computational  domain.  Notice  the 
impressive  size  and  strength  of  the  cells  on  the  right  hand 
side  of  the  storm.  At  times  their  long-axis  is  more  than 
100  m  with  pockets  of  downwelling  velocity  exceeding 
5u+  ~  0.35  ms1:  similar  large  values  were  observed  by 
vertical  profiling  floats  on  the  right  hand  side  of  Hurri¬ 
cane  Frances  (T.  B.  Sanford  2011,  personal  communica¬ 
tion).  Meanwhile,  on  the  left  hand  side  coherent  cells  can 
be  identified  but  they  are  clearly  weaker  in  magnitude 
and  smaller  in  scale.  This  flow  visualization  is  consistent 
with  the  temporal  evolution  of  vertical  velocity  variance 
and  third-order  moment  discussed  previously. 

Instantaneous  images  of  w  illustrate  the  chaotic  inter¬ 
action  between  cells  and  turbulence.  However,  in  order  to 
more  clearly  identify  the  average  coherent  structure  (a 
Langmuir  cell)  in  our  simulations  we  employ  the  technique 
of  linear  stochastic  estimation  (LSE)  (e.g.,  Adrian  et  al. 
1989;  Adrian  1996;  Christensen  and  Adrian  2001).  LSE  is 
a  robust  technique  for  identifying  structures  in  turbulent 
flows  and  is  an  excellent  approximation  of  a  conditional 
average  that  utilizes  only  unconditional  two-point  correla¬ 
tion  data.  For  our  flows,  we  estimate  the  deviation  from  the 
mean  horizontal  currents,  that  is,  the  horizontal  perturba¬ 
tion  velocities  ui  =  (w',i/,0)  at  x'  =  x  +  Ax  where  Ax  is 
the  spatial  displacement,  by  the  conditional  average 


(uV(x')lw(x)) 


(g(xK(x')) 

(Tv(x)w'(x)) 


w(x). 


(12) 


We  choose  resolved  vertical  velocity  W(x)  at  the  vertical 
location  of  (w2)max,  that  is,  z  ~  zmax,  as  the  event  indicator. 
LSE  utilizes  unconditional  two-point  spatial  correlations 
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FIG.  10.  Contours  of  vertical  velocity  normalized  by  u+  at  the  time 
of  maximum  winds  at  (top)  z  *  -7.2  m  along  the  track  X  -  -60  km 
and  (bottom)  at  z  —  — 12.6  m  along  the  track  X  =  60  km.  The 
subarea  depicted  in  the  top  (bottom)  panels  is  44  (11)  %  of  the  total 
x-y  area  in  their  respective  computational  domains.  The  black  (red) 
vectors  are  the  directions  of  the  surface  wind  (Stokes  drift). 


,  vgaa.  (i3a) 


^(ZmaxK  (*) 
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where  a  denotes  the  root-mean  square  of  the  given 
quantity.  In  our  implementation  of  (12)  and  (13)  with 
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FIG.  11.  Normalized  two-point  spatial  correlation  Rww(x\  y\  zmax)  along  the  track  X  -  60  km.  The  solid  (dotted)  contour  levels  are  0.1 
(-0.05).  The  positive  contour  level  encloses  a  region  of  strong  downwelling  while  the  negative  contour  level  encloses  a  broad  region  of 
weaker  upwelling.  The  vectors  (arrows)  show  the  magnitude  and  direction  of  the  horizontal  perturbation  velocities  («',  u')  induced  by 
a  Langmuir  cell,  (a)-(f)  t  -  tm  -  (-30.2,  -12.5,  -3.04,  -0.003, 3.37, 12.8)  h,  repectively.  Note  the  change  in  the  horizontal  and  vertical 
scales  between  the  lower  and  upper  panels. 

LES  data  the  angle  brackets  denote  an  average  over  an  with  the  findings  of  McWilliams  et  al.  (1997)  and  Sullivan 
x  —  y  plane  for  each  particular  point  (x',y')-  et  al.  (2007),  and  are  consistent  with  the  CL2  instability 

Time  sequences  of  the  conditionally  averaged  field  ui  mechanism  described  by  Leibovich  (1983).  Two  of  the 

at  z  =  zmax  and  contours  of  the  normalized  two-point  most  prominent  features  of  the  cells  in  the  hurricane  OBL 

spatial  correlation  Rww(x',y’,  zmax)  are  provided  in  Figs.  11  are  their  rotation  and  scale  expansion  and  contraction  for 

and  12.  We  define  the  spatial  extent  of  the  coherent  varying  time  and  track  position.  The  flow  patterns  exhibit 

Langmuir  structure  by  contours  of  Rww ;  the  positive  strong  horizontal  convergence  perpendicular  to  the  major 

contour  level  0.1  encompasses  an  ellipsoidal  region  of  axis  of  the  cell  center.  There  is  vigorous  inflow  at  the  cell 

strong  downwelling  with  side  lobes  of  broad  weaker  tail  and  weaker  outflow  at  the  cell  head  with  the  strongest 

upwelling,  contour  level  -0.05.  Inside  a  downwelling  perturbation  velocities  centered  just  above  the  largest 

region,  the  perturbation  velocity  field  ui  is  observed  to  downwelling  motions.  On  the  left-hand  side,  the  cells  ro¬ 
be  closely  aligned  (parallel)  to  the  major  axis  of  the  el-  tate  continuously  in  counterclockwise  fashion.  Close  in- 

liptical  shaped  structure.  The  coherent  structure  defined  spection  of  the  flow  pattern  at  t  —  tm  ~  0  shows  the 

by  the  two-point  correlation  function  Rww  and  the  ac-  development  of  a  vertical  rotational  motion  that  displaces 

companying  perturbation  velocity  fields  are  a  consistent  the  perturbation  velocity  vectors  slightly  off  center.  Van 

average  representation  of  the  instantaneous  flow  fields  Roekle  et  al.  (2012),  in  idealized  LES  with  misaligned 

shown  previously.  The  flow  patterns  are  in  good  agreement  winds  and  waves,  finds  a  similar  result  and  attributes  this 
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Fig.  12.  The  normalized  two-point  spatial  correlation  y',  zmax)  and  perturbation  horizontal  velocity  (m\  i/)  along  the  track 

X  =  -60  km:  contour  levels  and  labeling  as  in  Fig.  11;  (a)-(f)  correspond  to  t  —  tm  —  (  —30.9,  —12.7,  —3.44,  -0.003,  2.55,  10.6)  h, 
respectively.  Note  how  the  Langmuir  cells  rotate  counterclockwise  in  contrast  to  those  in  Fig.  11. 


to  advection  of  horizontal  vorticity  across  the  cell.  On  the 
right  hand  side  of  the  storm,  the  cells  first  rotate  counter 
clockwise  then  rapidly  clockwise.  The  cells  grow  rapidly  in 
scale  near  the  time  of  maximum  winds  as  they  tend  to 
align  themselves  with  the  surface  winds.  In  contrast,  the 
cell  growth  is  much  more  modest  on  the  left  hand  side  of 
the  storm,  consistent  with  the  temporal  changes  in  the 
vertical  velocity  moments.  On  the  left  hand  side,  the  cells 
are  larger  at  t  —  tm  =  (”12.7,  -3.44)  h  than  they  are  at  the 
time  of  maximum  winds. 

c.  Langmuir  cell  alignment 

One  of  the  important  questions  is  the  orientation  of 
a  Langmuir  cell  for  misaligned  winds  and  waves.  Cell 
orientation  significantly  complicates  the  development  of 
OBL  parameterizations  with  wave  effects  (McWilliams 
et  al.  2012).  In  Fig.  13,  we  show  the  surface  wind  and  Stokes 
drift  velocities  and  the  direction  of  the  Lagrangian  shear 
given  by  (11).  The  orientation  of  the  cells  in  LES  is  deduced 
from  the  conditionally  averaged  velocity  fields  sampled  at 
the  cell  center  [h'(0,  0,  Zmax),  i/(0, 0,  zmax)]*  The  Langmuir 


shear  Cs  is  a  key  quantity  in  Langmuir  turbulence  as  it  ap¬ 
pears  as  a  production  term  in  the  turbulent  kinetic  energy 
equation  [see  Eq.  (5.1)  from  McWilliams  et  al.  (1997)  and 
Grant  and  Belcher  (2009)].  In  addition,  Holm  (1996)  shows 
how  the  Lagrangian  shear  contributes  to  vortex  stretching, 
that  is,  the  tilting  of  vertical  vorticity  into  streamwise  di¬ 
rections.  The  results  in  Fig.  13  indicate  that  the  cells  are 
oriented  primarily  with  Cs ,  which  is  also  broadly  parallel  to 
the  surface  wind  direction.  The  deviation  from  £s  is  largest 
along  X  =  60  km  in  the  wake  where  the  cell  alignment 
tends  more  toward  the  wind  direction.  Figure  14  shows 
how  Cs  evolves  with  time  at  the  water  surface.  Notice 
that  there  are  periods  when  Stokes  and  current  shear  are 
of  opposite  sign,  which  leads  to  a  reduction  in  turbulence 
variances  (see  discussion  in  section  5a). 

6.  Impacts  of  Langmuir  turbulence  on  OBL 
statistics 

In  the  hurricane  driven  OBL,  coherent  Langmuir  cells 
interacting  with  background  turbulence  leave  an  imprint  on 
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FIG.  13.  Time  variation  of  the  surface  wind  (dotted  black  line),  the  surface  Stokes  drift  ve¬ 
locity  (gray  line),  and  the  Lagrangian  shear  (black  line):  (left)  X  -  -60  km  and  (right)  X  - 
60  km  .The  dominant  Langmuir  cell  direction  deduced  from  the  LES  are  indicated  by  solid 
circles.  Angles  are  measured  counterclockwise  from  the  x  axis. 


the  vertical  velocity.  Langmuir  turbulence  also  effects  the 
OBL  low-order  statistics,  namely,  mean  currents,  temper¬ 
ature  profiles,  momentum  and  heat  fluxes,  variances  and 
dissipation.  This  is  an  overly  abundant  amount  of  informa¬ 
tion  to  present,  and  thus  we  focus  on  a  few  select  statistics 
that  demonstrate  the  role  of  wave-current  interactions. 

a .  Horizontal  variances 

In  simulations  with  steady  aligned  winds  and  waves 
the  coherent  Langmuir  cells  modulate  the  horizontal 
streamwise  and  spanwise  variances  compared  to  their 
flat-wall  counterparts.  The  streamwise  oriented  cells  in¬ 
duce  strong  spanwise  convergence  at  downwelling  sites 
and  their  vortical  rotation  about  the  x  axis  amplifies  the 
spanwise  variance  (i/2)  at  the  expense  of  streamwise 
variance  (w'2).  For  example,  McWilliams  et  al.  (1997) 
reports  that,  based  on  the  resolved  LES  velocities,  the 
ratio  (u'2)/(m'2)  >3  near  the  water  surface;  theoretical 
linearized  rapid  distortion  results  for  a  propagating  sur¬ 
face  wave  are  in  qualitative  agreement  with  this  numer¬ 
ical  result  (Teixeira  and  Belcher  2002).  Wave  breaking, 
however,  tends  to  weaken  this  variance  anisotropy  by 
disrupting  the  near-surface  cell  structure  (Noh  et  al.  2004; 
Sullivan  et  al.  2007).  This  intercomponent  exchange  of 
horizontal  variance  is  also  present  in  our  simulations  and 
reflects  the  directional  changes  in  the  time  varying  winds 
and  waves  as  depicted  in  Fig.  15.  As  shown  previously, 
cells  tend  to  align  with  the  time  dependent  Lagrangian 
shear  direction  and  the  rotating  cells  induce  a  smooth 
exchange  of  variance  between  the  horizontal  velocity 


components.  The  cells  also  noticeably  alter  the  phase 
relationship  between  the  resolved  horizontal  variances 
so  that  (a'2)  and  (u'2)  reach  maximum  values  at  times 
markedly  different  than  simulations  without  vortex  force. 
This  is  especially  apparent  at  the  time  of  maximum  winds 
on  the  right  hand  side  of  the  storm  track.  In  the  absence  of 
vortex  force  (ut2)  is  a  minimum  at  t  —  tm  =  0,  while  it 
reaches  a  maximum  when  vortex  force  is  included.  Note 
that  at  t  -  tm  =  0  the  cells  are  primarily  oriented  in  the 


FIG.  14.  Mean  Lagrangian  shear  components  at  the  water  surface 
on  (top)  the  left-hand  side  and  (bottom)  right-hand  side  of  the 
storm.  The  shears  are  made  dimensionless  by  Zi/w*,  where  z\  is 
the  location  of  the  first  LES  gridpoint  below  the  water  surface;  the 
heavy  black  lines  indicate  time  periods  where  the  current  and 
Stokes  shears  have  opposite  signs. 
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Fig.  15.  Temporal  variation  of  the  resolved  horizontal  variances 
(h'2,U'2)  at  z  =  -0.5  m  below  the  water  surface  for  simulations 
with  and  without  wave  effects  at  locations  left  and  right  of  the 
storm  center:  simulations  with  wave  effects  for  (a)  X  =  60  km  and 
(b)  X  =  -60  km  and  with  no  wave  effects  for  (c)  X  =  60  km  and 
(d)  X  =  -60  km.  As  the  winds  and  waves  veer  with  time  the 
horizontal  component  with  maximum  variance  also  shifts  tracking 
the  rotation  of  the  coherent  Langmuir  cells.  At  t  -  tm  =  0  the  winds 
are  primarily  aligned  with  the  y  direction  and  then  the  v  component 
is  a  maximum  as  might  be  expected  for  a  wall-bounded  shear  flow 
without  wave  effects.  The  inclusion  of  vortex  force  alters  this  re¬ 
lationship  so  that  the  u  component  is  a  maximum. 

y  direction,  which  results  in  an  amplification  of  the  u 
variance  with  a  reduction  in  the  v  variance.  Our  flow 
visualization  shows  that  the  cells  are  weaker  and  smaller 
in  scale  on  the  left-hand  side  of  the  storm  track  and  this 
results  in  a  lower  value  of  u  variance  compared  to  its 
counterpart  on  the  right-hand  side  near  the  water  surface. 

b.  Current  and  variance  vertical  profiles 

Vertical  profiles  of  the  mean  currents  (Z7,u)  highlight 
the  complex  coupling  that  occurs  between  winds  and 
waves  and  between  waves  and  currents.  Figure  16  com¬ 
pares  the  mean  currents  and  Stokes  drift  velocity,  nor¬ 
malized  by  local  m*,  on  the  left  and  right  sides  of  the  storm 
track  for  simulations  at  the  time  of  maximum  winds. 
Wave  influences  are  pronounced.  On  the  right-hand  side 
the  Stokes  drift  velocity  is  dominated  by  the  vs  compo¬ 
nent  with  the  wind-Stokes  alignment  near  unity  over  a 
sizeable  fraction  of  the  OBL  (see  lower  panel  of  Fig.  5). 
The  potent  Langmuir  cells  induce  a  vertically  uniform 
u-current  profile  down  to  a  depth  zi\h\  ~  0.6.  In  contrast, 
in  the  simulation  without  wave  effects  the  vertical  gra¬ 
dient  d(u)/dz  is  surface  intensified  but  varies  continuously 
over  the  entire  OBL.  In  cases  with  and  without  wave 
effects,  there  are  steep  current  gradients  near  the  base  of 
the  OBL. 

Inspection  of  the  current  profiles  shows  that  Stokes 
drift  induces  a  return  Eulerian  flow  opposite  to  the 


direction  of  wave  propagation:  this  is  a  consequence  of 
Stokes-Coriolis  coupling, /z  X  nSj  in  the  LES  equations 
(Huang  1979;  Gnanadesikan  and  Weller  1995;  McWilliams 
et  al.  1997;  Polton  et  al.  2005).  On  the  right-hand  side  of 
the  storm  this  process  leads  to  Eulerian  current  changes 
that  are  similar  to  simulations  with  aligned  winds  and 
waves.  However,  on  the  left-hand  side  of  the  storm  with 
wind-wave  misalignment,  the  wave-induced  return  flow 
leads  to  oppositely  signed  gradients  d(u)ldz  and  dujdz . 
As  discussed  in  section  5a,  the  vertical  momentum  flux 
(u'w)  is  then  countergradient  to  the  gradient  of  Stokes 
drift,  which  inhibits  the  production  of  w  variance.  Over 
the  bulk  of  the  OBL  the  currents  with  and  without  wave 
effects  are  broadly  similar  with  the  important  changes 
occurring  near  the  water  surface. 

The  resolved  velocity  variance  profiles  displayed  in 
Fig.  17  paint  a  consistent  picture  of  the  impact  of 
Langmuir  turbulence  in  the  nonequilibrium  hurricane 
regime.  On  the  left-hand  side  of  the  storm  at  t  -  tm  =  0, 
the  w  variance  in  the  simulation  with  vortex  force  dis¬ 
plays  a  slightly  elevated  subsurface  maximum  compared 
to  its  counterpart  without  vortex  force.  The  wave  effects 
on  the  horizontal  variances  are  confined  to  a  shallow 
layer  near  the  surface  -0.15  <  zi\h\  <  0,  and  outside  this 
region  the  profiles  from  the  simulations  with  and  without 
vortex  force  are  nearly  indistinguishable.  This  is  consis¬ 
tent  with  the  flow  visualization  and  the  interpretation  of 
weak  Langmuir  turbulence  on  the  left  hand  side  of  the 
storm  discussed  in  section  5.  The  variance  profiles  along 
X  =  60  km  are  in  sharp  contrast  to  their  counterparts 
along  X  =  —  60  km.  The  wave  effects  are  pronounced  and 
the  changes  in  the  v  and  w  variances  extend  over  a  large 
fraction  of  the  OBL,  more  than  50%.  Also,  it  is  apparent 
that  vertical  transport  by  the  Langmuir  cells  contributes 
to  the  elevated  v  variance  in  the  region  -1  <  zi\h\  <  -0.5. 
This  hints  at  a  coupling  with  shear  instabilities  at  the 
thermocline  induced  by  inertial  resonance. 

c.  Entrainment  flux  and  SST 

One  of  the  primary  bulk  functions  of  the  hurricane 
OBL  is  to  entrain  cool  water  from  below  the  thermo¬ 
cline  into  the  mixed  layer.  Over  time,  this  lowers  SST 
that  can,  depending  on  the  translation  speed  of  the 
storm,  feedback  onto  the  large-scale  hurricane  dynam¬ 
ics.  Previous  modeling  and  observational  studies  (e.g., 
Price  1981;  Crawford  and  Large  1996;  Skyllingstad  et  al. 
2000;  Sanford  et  al.  2011)  find  that  the  primary  dynam¬ 
ical  mechanism  for  entraining  cool  water  is  shear  in¬ 
stability  at  the  base  of  the  OBL  induced  by  inertially 
rotating  wind  stress  (see  the  upper  panel  of  Fig.  5).  This 
critical  dynamical  process  is  also  found  to  be  dominant 
in  our  simulations  but  is  augmented  by  additional  mix¬ 
ing  induced  by  the  coherent  Langmuir  cells  that  develop 
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FlG.  16.  Profiles  of  mean  currents  (u,v)  and  Stokes  drift  velocity  (uSy  vs )  normalized  by  «*  at 
the  time  of  maximum  winds:  left  (right)  sides  of  the  storm  are  shown  in  left  (right)  panels.  The 
black  (gray)  colors  are  simulations  with  (without)  vortex  force.  Averages  are  computed  over  an 
x-y  plane  and  3600  s  in  time. 


in  the  OBL  on  both  sides  of  the  storm  track.  This  is  il¬ 
lustrated  in  Fig.  18  where  we  compare  the  minimum  total 
heat  flux  at  the  base  of  the  OBL,  that  is,  the  entrainment 
flux  (w>0')min,  and  the  change  in  SST,  A0  =  (0)  -  0,  at 
z  =  0,  for  simulations  with  and  without  wave  effects.  The 
entrainment  flux  (w>0')mjn  is  the  sum  of  LES  resolved  and 
SGS  contributions  and  is  simply  computed  by  sweeping 
through  the  vertical  profile  of  the  mean  total  temperature 
(or  heat)  flux  and  locating  its  minimum  value.  Vertical 
profiles  of  mean  temperature  and  heat  flux  are  obtained 
by  averaging  over  an  x  —  y  plane  in  the  LES  domain  and 
1000  s  in  time. 

Far  ahead  of  the  storm  center  the  normalized  entrain¬ 
ment  flux  (w&f)m\nIQ+(X,  i)  - - 0.25,  which  is  slightly 

more  negative  than  the  classical  value  of  -0.2  for 
a  daytime  convective  atmospheric  boundary  layer  (e.g., 
Schmidt  and  Schumann  1989;  Moeng  and  Sullivan  1994). 
This  is  expected  since  the  Stokes  drift  is  modest  and 
the  OBL  is  in  a  regime  of  mixed  shear  and  convective 

forcing  with  the  ratio  |/i,|/L - 1.  With  increasing  time, 

the  OBL  becomes  increasingly  dominated  by  shear  and 
wave  forcing.  At  the  time  of  maximum  winds,  the  en¬ 
trainment  flux  exceeds  its  surface  value  by  a  factor  of  (—10, 
-3)  on  the  (right,  left)  sides  of  the  storm,  respectively. 
Langmuir  turbulence  noticeably  enhances  the  entrain¬ 
ment  flux  compared  to  the  situation  without  wave 


effects.  For  example,  on  the  right-hand  side  at  t  —  tm  = 
-4  h  the  normalized  entrainment  flux  increases  from 
-4  to  -6,  a  50%  increase  when  wave  effects  are  in¬ 
cluded.  At  t  —  tm  =  0,  the  increase  is  about  20%  and  the 
resulting  maximum  change  in  A0  is  about  0.25  K  cooler. 
In  absolute  values,  the  changes  in  entrainment  flux  and 
SST  on  the  left  hand  side  of  the  storm  are  more  modest 
than  their  counterparts  on  the  right  hand  side.  On  both 
sides  of  the  storm  track  wave  effects  cause  the  minimum 
entrainment  flux  to  occur  earlier  in  time,  perhaps  by  as 
much  as  3  h.  This  suggests  enhanced  coupling  with  the 
inertial  currents  at  the  thermocline.  Our  estimate  that 
the  normalized  entrainment  (cooling)  flux  on  the  right- 
hand  side  of  the  storm  is  greater  than  10  times  the  sur¬ 
face  heat  flux  is  in  good  agreement  with  the  observations 
by  Sanford  et  al.  (2011,  see  p.  1044).  The  LES  predicted 
changes  in  SST  on  the  left  and  right  hand  sides  of  the 
storm  are  in  good  agreement  with  the  observations  of 
D’Asaro  et  al.  (2007),  but  those  observations  do  not 
readily  lend  themselves  to  quantifying  the  importance  of 
surface  waves. 

In  the  hurricane  wake  the  entrainment  flux,  especially 
on  the  right-hand  side  of  the  storm,  appears  to  retreat 
toward  its  prestorm  value  in  a  series  of  steep  rises  and 
declines  that  is  noticeably  different  than  in  the  time 
period  t  -  tm  <  0.  Flow  visualization  reveals  packets  of 
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FlG.  17.  (left  to  right)  Vertical  profiles  of  resolved  variances  {h'2,D'2,  w'2)  normalized  by  m*  for  (top)  left  side  and 
(bottom)  right  side  of  the  storm  track.  The  black  (gray)  solid  lines  in  each  figure  are  simulations  with  (without)  vortex 
force.  Results  are  at  t  -  tm  -  0. 


large-scale  propagating  gravity  waves  that  overturn  in¬ 
termittently,  creating  large  oscillations  in  the  entrain¬ 
ment  flux  as  well  as  other  OBL  statistics.  These  turbulent 
gravity  waves  appear  to  be  a  product  of  the  vigorous 
resonant  wind  stress  forcing.  The  large  horizontal  scale 
of  these  structures  motivated  our  use  of  the  wide  LES 
horizontal  domain  described  in  section  3b. 

d.  Boundary  layer  depth  and  Stokes  scale 

It  is  interesting  to  notice  that  the  augmentation  of  the 
bulk  temperature  entrainment  and  SST  change  induced 
by  Langmuir  turbulence  occurs  during  the  time  period 
t  —  tm<  0.  This  can  be  interpreted  in  terms  of  the  surface 
La,  and  the  variations  of  the  bulk  OBL  length  scales  h 
and  Ds .  Figure  19  depicts  the  time  growth  of  the  OBL 
and  the  ratio  of  the  Stokes  depth  scale  to  the  local  OBL 
depth,  that  is,  the  ratio  DJ\h\,  along  the  tracks  X  = 
±60  km.  The  OBL  depth  h  is  estimated  using  the  max¬ 
imum  temperature  gradient  method  described  in  Sullivan 
et  al.  (1998);  essentially  h  is  the  spatial  average  of  the 
boundary  layer  depths  /z(x,  y)  deduced  along  every 


vertical  column  in  the  LES  grid  box.  A  comparison  of 
Fig.  18  with  Fig.  19  shows  that  the  largest  impact  of 
Langmuir  turbulence  on  the  entrainment  flux  occurs  in 
the  deepening  phase  of  the  OBL,  that  is,  when  the  ratio 
DJ\h\  >  0.75;  a  similar  result  is  reported  by  Skyllingstad 
et  al.  (2000).  Based  on  the  definition  of  Ds  given  by  (7), 
this  empirical  ratio  corresponds  to  a  nondimensional  e- 
folding  length  scale  for  a  monochromatic  Stokes  wave  of 
zd\h\  =  (3/16 ?r)  or  about  6%  of  the  OBL  depth.  We 
note  that  in  the  simulations  of  McWilliams  et  al.  (1997) 
and  Sullivan  et  al.  (2007)  DJ\h\  ~  1.82  and  they  find 
significant  wave  effects  on  entrainment.  In  the  present 
simulations  when  DJ\h\  >  0.75,  the  coherent  cells  are  of 
sufficient  strength  to  effectively  transport  momentum 
and  kinetic  energy  to  the  base  of  the  thermocline  and 
thereby  increase  the  entrainment  of  cool  water.  In  the 
hurricane  wake,  t  —  tm  >  0,  the  ratio  DJ\h\  <  0.5  because 
the  OBL  mixes  to  very  deep  levels  by  inertial  resonance 
and  Ds  is  reduced  by  the  falling  surface  winds.  The  ratio 
DJ\h\  appears  to  be  an  empirical  indicator  of  when  Lang¬ 
muir  turbulence  impacts  entrainment  in  a  stratified  OBL. 
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Fig.  18.  Temporal  variation  of  (top)  entrainment  flux  and  (bottom)  SST  change.  The  groups  of 
curves  marked  “L”  (“R”)  are  conditions  on  the  left  (right)  hand  side  of  the  storm.  In  each  figure 
the  black  (gray)  lines  are  simulations  with  (without)  vortex  force.  The  entrainment  flux  is  nor¬ 
malized  by  Q*(X,  t)  the  local  value  of  the  surface  heat  flux  (see  Fig.  4)  and  A 0  =  (0)  -  0,-  at  z  -  0. 


7.  Summary  and  conclusions 

A  turbulence  resolving  large-eddy  simulation  (LES) 
model  is  used  to  study  ocean  boundary  layers  (OBLs) 
driven  by  extreme  (hurricane)  winds  and  waves.  The  LES 
model  adopts  wave-averaged  equations  that  include  the 
vortex  force  us  X  to,  where  uv  is  the  wave  field  Stokes  drift 
velocity  and  a>  is  the  vorticity.  Large-scale  time-varying 
wind  and  wave  fields  input  to  the  LES  are  generated  by 
WaveWatch  III;  the  wind  and  wave  fields  are  an  ide¬ 
alization  of  Hurricane  Frances  traveling  at  a  speed  of 
5.5  m  s“\  The  surface  wind  and  wave  fields  are  in  a  non¬ 
equilibrium  state  and,  as  a  result,  the  Stokes  drift  veloc¬ 
ity  is  a  complex  two-component  vector  that  varies  with 
time,  storm  track  position,  and  spirals  with  depth.  High- 
resolution  LES  using  as  many  as  1024  X  1024  X  256  grid 
points  integrated  for  more  than  70  physical  hours  are 
conducted  at  fixed  positions  right  and  left  (X  =  ±60  km) 
of  the  storm  track. 

The  surface-based  turbulent  Langmuir  number  La„ 
which  is  used  to  assess  the  relative  importance  of  shear 
and  wave  forcing,  varies  from  0.2  to  more  than  0.6  but 
over  wide  regions  La,  —  0.3.  Thus  the  hurricane  forced 
OBL  is  often  in  a  Langmuir  turbulence  regime.  The  LES 
solutions  exhibit  important  differences  from  previous 
simulations  that  adopt  wind  wave  equilibrium.  The  po¬ 
tency  of  Langmuir  turbulence  and  the  size  and  strength 


of  the  coherent  Langmuir  cells  depends  on  the  locally 
evolving  wind  wave  state.  At  the  time  of  maximum  winds 
on  the  right-hand  side  of  the  storm,  the  near  inertial 
rotation  of  the  winds  creates  a  “resonant”  Stokes  drift 
that  is  well  aligned  with  the  surface  winds  over  a  large 
fraction  of  the  OBL.  Imagery  shows  that  the  resulting 
instantaneous  Langmuir  cells  can  be  more  than  100  m  in 


FIG.  19.  Temporal  variation  of  the  OBL  depth  \h\  (left  ordinate) 
and  ratio  DJ\h\  (right  ordinate).  The  solid  black  (gray)  lines  are 
simulations  with  vortex  force  on  the  right  (left)  side  of  the  storm. 
The  Stokes  depth  scale  Ds  is  predicted  by  (7).  The  dashed  black 
(gray)  lines  are  estimates  of  \h\  on  the  right  (left)  side  for  simula¬ 
tions  with  no  vortex  force.  For  reference,  the  LES  by  McWilliams 
et  al.  (1997)  is  carried  out  with  the  ratio  DJ\h\  ~  1.82. 
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length  with  downwelling  velocities  exceeding  5m*  (m*  is 
the  friction  velocity)  and  are  depth  filling.  The  cells  are 
much  weaker  and  smaller  on  the  left-hand  side  because 
of  wind  wave  misalignment.  The  latter  results  in  oppo¬ 
sitely  signed  vertical  gradients  of  current  and  Stokes  drift. 
Hence  the  vertical  momentum  fluxes  can  be  counter- 
gradient  to  dnjdz ,  which  results  in  less  production  of 
vertical  velocity  variance  near  the  water  surface.  Verti¬ 
cal  momentum  fluxes  countergradient  to  vertical  gradi¬ 
ents  of  Stokes  drift  are  one  of  many  possible  regimes  in 
nonequilibrium  winds  and  waves.  The  average  coherent 
Langmuir  cell,  identified  using  linear  stochastic  estima¬ 
tion,  rotates,  grows,  and  shrinks  depending  on  La„  and  the 
alignment  between  the  surface  wind  and  depth-varying 
Stokes  profile.  Over  time  on  both  sides  of  the  storm  track, 
the  cells  tend  to  closely  align  with  the  mean  Lagrangian 
shear  direction. 

The  vertical  profiles  of  mean  currents  and  tempera¬ 
ture,  turbulence  variances,  and  momentum  and  heat 
fluxes,  are  impacted  by  Langmuir  turbulence  to  varying 
degrees.  At  the  time  of  maximum  winds,  Langmuir  tur¬ 
bulence  enhances  the  minimum  heat  flux  (or  entrainment 
flux)  at  the  base  of  the  OBL  by  20%  and  lowers  the  sea 
surface  temperature  by  0.25  K  on  the  right-hand  side  of 
the  storm.  Inertially  rotating  wind  stress,  however,  re¬ 
mains  the  dominant  dynamical  process  for  entraining 
cold  water  into  the  OBL.  Wave  effects  also  produce 
changes  to  entrainment  and  SST  on  the  left-hand  side 
of  the  storm,  but  to  a  lesser  degree.  In  the  non¬ 
equilibrium  wind  wave  regime  studied  here  the  surface 
based  turbulent  Langmuir  number  La,  appears  to  be 
insufficient  in  describing  when  wave  effects  impact  the 
OBL.  The  shape  of  the  Stokes  drift  vertical  profile, 
wind  wave  alignment,  and  the  depth  scale  of  the  Stokes 
profile  Ds  are  important  bulk  parameters.  When  the 
ratio  DJ\h\  >  0.75,  where  h  is  the  OBL  depth,  Langmuir 
turbulence  enhances  entrainment  at  the  thermocline. 
Wave  breaking,  which  is  a  dominant  feature  of  the  high 
wind  OBL,  is  omitted  here  but  will  be  considered  in 
future  work. 
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